Cosmic string loops and large-scale structure 
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We investigate the contribution made by small loops from a cosmic string network as seeds for 
large-scale structure formation. We show that cosmic string loops are highly correlated with the 
long-string network on large scales and therefore contribute significantly to the power spectrum of 
density perturbations if the average loop lifetime is comparable to or above one Hubble time. This 
effect further improves the large-scale bias problem previously identified in earlier studies of cosmic 
string models. 



A. Introduction 

Quantitative predictions for the large-scale structure 
induced by cosmic strings have taken some time to crys- 
tallise as the understanding of cosmic string physics has 
improved [El. In particular, the role of small loops pro- 
duced by the string network has evolved from a potential 
one-to-one correspondence between loops and cosmolog- 
ical objects Q through to a completely subsidiary role 
relative to the wakes swept out by long strings ||] . This 
dethronement of loops was a result of numerical stud- 
ies which showed that the average loop size I = at was 
much smaller than the horizon, I « g?h (§,[§; ^hcy might 
even be as small as the lengthscale set by gravitational 
backreaction a. ~ 10 -4 , a value appropriate for GUT- 
scale strings M. Add the high ballistic loop velocities 
observed v ss c/\/2 and it was not surprising that these 
tiny loops have been assumed to be more or less uni- 
formly distributed and hence a negligible source relative 
to the long string network ||. Nevertheless, small loops 
always make up a significant fraction of the total string 
energy density at any one time and, as we demonstrate 
here, loop-induced inhomogeneities are considerable if 
their lifetime is not much smaller than the Hubble time. 
By properly incorporating these loop perturbations, we 
show that their contribution relative to the long string 
wakes is almost comparable and also highly correlated 
with these wakes. 

The context for this work is a major programme of 
structure formation simulations seeded by high resolution 
cosmic string networks with very large dynamic ranges 
. This work demonstrated that for open or A mod- 
els with r = Vlh — 0.1-0.2 and a cold dark matter 
(CDM) background, the linear density fluctuation power 
spectrum has both an amplitude at 8/i _1 Mpc, as, and 
an overall shape which are consistent within uncertain- 
ties with those currently inferred from galaxy surveys. 
This result has also been confirmed using semi-analytical 
phenomenological models which incorporated some of the 
main features of long string networks JlO| , |ll[ . 

In this letter we investigate the contribution of cosmic 
string loops to the linear power spectrum of cosmic string 



induced density perturbations. This component has been 
ignored and excluded in previous work, due to both the 
computational difficulties and assumptions about the ho- 
mogeneity of the loop distribution. To this end we first 
perform very high resolution numerical simulations of a 
cosmic string network with a dynamic range extending 
from well before the radiation-matter transition through 
to deep into the matter era. We then use this net- 
work as a source for density perturbations (as described 
in taking into account the large-scale power con- 

tributed from cosmic string loops. This is done by mod- 
elling cosmic string loops smaller than a fixed fraction of 
the horizon size as relativistic point masses. The effects 
of the evaporation of these loops into gravitational waves 
and the damping of loop motion due to expansion are 
also included. Note, however, that this should be clearly 
distinguished from recent work fl^| , which attempts to 
incorporate network decay products in the power spec- 
trum of an additional fluid (with a variety of possible 
equations of state). This does not appear to properly ac- 
count for the phase correlation between long strings and 
moving loops. 

Unless otherwise indicated, we use h — 0.7, £l m = 1 
and f^A = in the results presented here. A verified ac- 
curate rescaling scheme for the resulting power spectrum 
with different choices of h, f2 m and J7a is straightforward 
and described in ref . |7],||,[l(^[l3| . 

B. Cosmic string and loop evolution 

The Nambu equations of motion for cosmic strings in an 
expanding universe can be averaged to yield: 

^ + 2H(l + (v 2 )) Poo = -X L , (1) 

where is the long string energy, t the physical time, 
H = a/ a the Hubble parameter, a(t) the scale factor, (v 2 ) 
the mean square velocity of strings, and the transfer 
rate of energy density from long strings into loops. In 
the scaling regime the long string energy density should 
scale with the background energy density evolving as 
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dpoc 2 P°° 

dt ~ T' 

Substituting this into (0) to eliminate dpoo/dt gives 
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(1 — (w 2 )) ~ 0-6 in radiation era, 
|(1 - 2(?4)) - 0.2 in matter era, 



(2) 



(3) 



where (v 2 ) > «) - 0.6 §,§. Both (|) and (|) provide 
a check for the scaling behavior of long strings and loops 
in the cosmic string network simulations. 

We know that the loops produced by a cosmic string 
network will decay into gravitational radiation, with a 
roughly constant decay rate TGp 2 , where fi is the string 
linear energy density. Typically T = 50— 100 with an av- 
erage (r) ~ 65 ]li|Jl^| . Now, if we assume the loop pro- 
duction to be 'monochromatic' so that all loops formed 
at the same time will have the same mass, we can write 
the initial rest mass of a loop formed at time as 



Ml = apt* = fTGfiH, 



(4) 



Here, the parameter / = a/TGfi is expected to be of or- 
der unity if the size of the loops formed at the time t is de- 
termined by gravitational radiation back-reaction, which 
smoothes strings on scales smaller than TG/it. With the 
decay rate introduced earlier, we have the rest mass of a 
loop formed at time i* evolving as 



M L {U,t) =MlW(U,t) , 



(5) 



where 
W{U,t) 



t-t, 
r(n) 

otherwise 



1 - ^jy for t*<t<U + t{U) 



Here r(t^) « /t* is the lifetime of loops produced at 
time t* (/ = 2, 3 implies the decay occurs in one horizon 
time in the radiation and matter eras respectively). The 
evolution of the loop energy density is then given by: 



PL(t) = 




ait') 



1 3 



W(t',t)dt' 



a(t) 

for / < 1, 

for / ^> 1 (radiation era), 
(ln/)/i 2 for / 3> 1 (matter era), 



(7) 



where we have used the scaling behaviour (||) and ||). 
Consequently, the scaling of the power spectrum induced 
by loops in / should interpolate between f 2 and / (radi- 
ation era) or (In/) 2 (matter era). We notice in (Q) that 
we have ignored the effect of loop velocity redshifting due 
to the expansion of the Universe, which causes a change 
in the effective mass. Because loops are formed with rel- 
ativistic velocities, we expect this damping mechanism to 
have the strongest effect for / 3> 1, but to be negligible 
for / < 1. 
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FIG. 1. Evolution of X L (t). The dotted lines are the 
asymptotic values in the radiation and the matter eras. 
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If a loop formed at time t* has an initial physical ve- 
locity v*, its trajectory in physical space accounting for 
the expansion of the Universe is then given by: 



x(t) = x.(U) +a(t) 



dt' 



t, y/a(t') 2 + A 2 a(t') 



(8) 



for t > t*, where A = 7*v*a*, A — |A| and 7* = 
(1 — |v» I 2 ) 1 / 2 . Here we have neglected the acceleration of 
loops due to the momentum carried away by the gravita- 
tional radiation, the so-called 'rocket effect'. A numerical 
calculation for several asymmetric loops shows that the 
rate of momentum radiation from an oscillating loop is 



TpGfi 2 



(9) 



where Tp ~ 10 |L6| |. Combining with (|^), one can show 
that this rocket effect will become important only when: 



t 



> 1 



1 + !>/(!» 
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which affects only the final stages of the loop lifetime as 
long as Tp/(T^v) < 1, or equivalently v > 0.15c. For a 
typical ~ c/%/2, one requires a loop lifetime > 43£* 
in the radiation era and > 17i* in the matter era to 
redshift down to this critical velocity according to (||). 
Since the values of / we explore here are of order unity, 
it is a reasonable approximation to neglect the transfer 
of momentum due to gravitational radiation. 

C. Results and discussion 

We first perform string simulations with a string sam- 
pling spacing 1/1000 of the simulation box sizes. The 
dynamic ranges cover from 0.05 to 300 rj c(l , where r] eq 
is the conformal time at radiation-matter energy density 
equality. We then perform the structure formation simu- 
lations with box sizes ranging from 20-120/i _1 Mpc, and a 
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FIG. 2. Small dynamic range power spectra of density per- 
turbations seeded by long strings (thick solid), by loops with 
initial velocities v, switched to zero (dot-dashed), and by 
loops with v* determined by string network evolution (dashed 
and thin solid). The thin solid line includes the effect of grav- 
itational decay of the loop energy, while the other two loop 
lines don't but with loops removed after a period of time 
r» = ir- 
resolution of 128 3 -512 3 . Figure [j] shows the evolution of 
Xl- We can see that the expected amount of energy was 
converted into loops in our simulations so that Xl has the 
correct asymptotic behavior given by ([5j). However, the 
typical loop-size (and consequently their lifetime) does 
not approach scaling so rapidly and is therefore larger 
than physically expected for most of the duration in the 
simulations To overcome this problem we rescale 

the loop lifetime according to equation (^). Thus the un- 
certainty in the average mass and therefore the lifetime 
of loops formed at a given time is quantified by the choice 
of the parameter /. The initial rms velocity of loops ob- 
served from the simulations is (v 2 ) 1 / 2 > 0.7c throughout 
all the regimes. 

Figure [2| shows the power spectrum of density pertur- 
bations induced by long strings and by cosmic string 
loops for / = 1 for a small dynamic range from 2.5 to 
5f? q- We can see that when compared with the spec- 
trum induced by static loops (dot dashed) , the amplitude 
of small-scale perturbations induced by moving loops 
(dashed) is clearly reduced by their motion. However, 
their large-scale power is higher because of the depen- 
dence of the gravitational interaction on the loop ve- 
locities, especially when they are relativistic. We also 
see that the gravitational decay of loop energy (thin 
solid) damps the overall amplitude of the power spec- 
trum (dashed) by about a factor of 3. We notice that 
between the long-string correlation scale ss 2O/77 
and the scale &l ~ 10/cj, the slope of the loop spec- 
trum (thin solid) is exactly the same as that of the long- 
string spectrum n sa —2.25 Q. We believe that this 
close correspondence is due to copious loop production 
being strongly correlated with long string intercommut- 
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FIG. 3. The lower set of 3 lines are for / = 0.5 (dot- 

ted), 1 (dot-dashed) and 2 (dashed). Voo(k) is plotted as a 
solid line. The upper set of lines are Vtot(k) with correspond- 
ing line styles and / values to the lower set of lines. 
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FIG. 4. The correlation coefficient between the long-string 
and loop induced perturbations, with / = 0.5, 1, 2, 4, 6 (down- 
wards) . 

ing events and the collapse of highly curved long string 
regions j|, that is, near the strongest long string per- 
turbations. Moreover, these correlations persist in time 
with the subsequent motion of loops and long strings ly- 
ing preferentially in the same directions, a phenomenon 
which has been verified by observing animations of string 
network evolution. These correlations between loops and 
long strings, however, have a lower cutoff represented by 
the mean loop spacing ~ fc£~ . Below d^, the effects 
of individual filaments swept out by moving loops can be 
identified. In terms of the power spectrum, for k < 
the loops are strongly correlated with the long strings and 
therefore reinforce the wake-like perturbations, while for 
k > /cl their filamentary perturbations increase the spec- 
tral index by about one to n w —1.25; this change is 
expected on geometrical grounds. 

In figure || we plot the power spectra of density per- 
turbations seeded by long strings "Poo(^), by small loops 
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FIG. 5. Comparison of the observational power spectrum 
@ with Voo{k) (solid), and P to t(fc) for / = 0.5,1,2,4,6 
(dashed, upwards), with a full dynamic range. 

Vh(k), and by both loops and long strings Vtot(k). The 
dynamic range here extends from 0.6 to 7.5?7 C q. As ex- 
pected Vh(k) scales more moderately than / 2 but more 
strongly than / (see (0)). It is also apparent that the per- 
turbations induced by long strings and by loops are posi- 
tively correlated with Vtot(k) > "Ph(k) + 'P ocl (k) through- 
out the whole scale range. This positive correlation 
between loops and long strings boosts the large-scale 
Voo{k) by a factor of 1.5, 1.8 and 2.2 to reach Ptot{k) 
for / = 0.5, 1 and 2 respectively, even if Vh{k) is a rela- 
tively small fraction of Voo(k) on these scales. 

Figure [| shows the correlation coefficient between the 
long-string and loop induced perturbations. We see that 
long strings and loops are strongly positively correlated 
on large scales, but weakly correlated on small scales 
where the loops dominate the perturbations (also see fig- 
ure ||) . The threshold k t between these two regimes must 
be significantly larger than /c L because, for k < k L , Ph(k) 
is well below and roughly parallel to Voo(k) (see figures 
^ and ||). We also verify that Vtat{k) /V^lk) is approxi- 
mately a constant for k < k^ <k t , which again provides 
strong evidence for the fact that loops behave as part of 
the long-string network on large scales. 

Given these properties of the string power spectra, one 
can easily construct a semi-analytic model for 'Ptot(fc) as 
for Voc{k) . We first multiply the structure function 
T(k,r]) of Poo{k) by J{f],f) to account for the boost 
'Ptot(k) /Poo(k) on large scales (k < fc t ), and then by a 
numerically verified form Ti(k, r], f) = [1 + (k/kt) 4 ] 1 ^ 4 to 
account for the turnover for k > k t (rj, /). J{r\, /) is cali- 
brated phcnomcnologically from simulations deep in the 
radiation era through to those deep in the matter era. In 
the pseudo-scaling regime for the loop size, kt is revealed 
to be at least 10fc^ 200/ry depending on /. Thus we 
can carry out a full-dynamic-range integration to obtain 
Vtot(k)- In figure || we compare this Vtot(k) and V x (k) 
j^U with observations |l7j]. The background cosmology 
is il c — 0.15, fl\ — 0.85 and h = 0.7, and we have used 



the COBE normalization Gfi = 1.7 x 10~ 6 JT|] through- 
out. Since loops are point-like and they have little im- 
pact through the Kaiser-Stebbins effect on COBE-scale 
CMB anisotropics, we expect this normalization to be 
very weakly dependent on the value of /; indeed, loops 
were found to be negligible in ref. |llj . Thus we see from 
figure | that for / > 0.5, loops can contribute signifi- 
cantly to the total power spectrum and ease the large- 
scale bias problem seen previously |],|9|,^0| . Definite con- 
clusions, therefore, about biasing in cosmic string models 
will need further advances in determining the magnitude 
of the parameter /, while all future large-scale structure 
simulations will now require the the inclusion of loops. 

These additional complications in modelling cosmic 
string structure formation are most obvious on small 
scales, where even higher resolution and large dynamic 
range simulations will be required. However, we expect 
the power spectrum on large scales to be only weakly 
dependent on the details of loop formation. Within the 
present pseudo-scaling regime for loop size, we know that 
kt > > as shown in figure ^| and discussed pre- 
viously. Taking this extreme minimum fc t = 10fc^, then, 
we find that the semi-analytic model over the full dy- 
namic range gives at most a 2% difference in Vtotik) 
for k < 1/iMpc -1 when the filament term Tt(k,r],f) 
is excluded from T(k,ri) (for il c = 0.15, Qa = 0.85 
and h = 0.7). This means that although the simula- 
tions described in this letter are already on the verge of 
present computer capabilities, a further detailed study on 
small scales will improve only the overall normalization 
of Vtot(k) but not the shape revealed here, which should 
be a robust feature. We note that advances in under- 
standing loop formation mechanisms will also be crucial 
in quantifying the importance of the gravitational radi- 
ation background emitted by cosmic string network and 
its effect on large-scale structure and CMB anisotropics 




D. Conclusion 

In this Letter we have described the results of high- 
resolution numerical simulations of structure formation 
seeded by a cosmic string network with a large dynamic 
range, taking into account for the first time the loops 
produced by the network. We show that on large scales 
the loops behave like part of the long-string network and 
can therefore contribute significantly to the total power 
spectrum of density perturbations, provided their life- 
time is not much smaller than one Hubble time. At 
present, the typical size and lifetime of loops formed by 
a string network remains to be studied in more detail; 
the problem is both computationally and analytically 
challenging. However, within the scale range of inter- 
est further developments in this area have the potential 
to affect the overall amplitude of the spectrum, while 
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leaving the shape largely unchanged. The results pre- 
sented here provide further encouragement for more de- 
tailed work on both the nature of cosmic string evolution 
and the large-scale structures they induce in cosmologies 
with r = Qh = 0.1-0.2. 
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